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. . . We study dynamics of electrons in a magnetic field using a network model with two channels 

^ : per link with random mixing in a random intrachannel potential; the channels represent either two 

Landau levels or two spin states. We consider channel mixing as function of the energy separation of 
the two extended states and show that its effect changes from repulsion to attraction as the energy 
separation increases. For two Landau levels this leads to level floating at low magnetic fields while 
for Zeeman split spin states we predict level attraction at high magnetic fields, accounting for ESR 
data. We also study random mixing of two degenerate channels, while the intrachannel potential is 
periodic (non-random). We find a single extended state with a localization exponent ~ f .f for real 
' scattering at nodes; the general case has also a single extended state, though the localized nature 

of nearby states sets in at unusually large scales. 
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^ ! I. INTRODUCTION 

1^ . 

The Quantum Hall Effect (QHE) remains a great attraction for theoreticians and experimentalists. Of particular 
' interest is the divergence of the localization length at a discrete set of energies, corresponding to extended states. 

a This has motivated a large variety of methods for studyingiihe pertinent metal insulator transition as a cri|tjcal second 
I ' order phase transitionEJ, such as field theoretical methodsEI, semiclassical methodscra, numerical methoda3 and finite 
[ size scaling for transfer matricesQ. Most of these works focus on properties near an isolated extended state with 
^ ■ emphasis on the critical exponent v of the Jocalization length. The most recent experimental resultQ v = 2.3 ±0.1 is 
O [ in agreement with theoretical predictionsO'E . 

, The situation with a few extended states, allowing for coupling between these states received less attention. This 

situation is relevant to the behavior of delocalized states in weak magnetic fields, ft is well known that the existence 
^— I , of delocalized states is a necessary condition for the QHE behavior. On the other hand, scaling theoryu and numerical 
K>" ' studiesB imply that in a 2D system in the absence of a magnetic field all states should be localized. Therefore , 
a scenario-iH-| which delocalized states "fioat up" above the Fermi level as the magnetic field decreases has been 
' suggestedMtil. Recent experimentala show that, indeed, the energy of of the lowest delocalized state floats up above 
[ the Fermi level as magnetic field is reduced. This corresponds to a transition frcagjaji insulator to a quantum Hall 
I conductor at both low and high fields. Numerical studies of a few Landau bandal3~llZl focused on critical exponents 
\^ ' except for an early work by Anddj which supports the floating scenario. Ando's work used S function impurities 
^\ , which is, however, not suitable when the impurity_concentration is too low since bound states of the S potential shift 
the extended state even for a single Landau bandO. We have recently shownE^I that two extended states attract each 
^ other which leads to a miniHium in the energy of the lower state. This attraction has been also predicted theoretically 
by Shahbazyan and Raikhllj using a high magnetic fleld expansion. _ _ 

[ The system of two coupled extended states was also studied in the context of a spin-split Landau bandHrtZl. In the 
absence of a Zeeman term it was found that two separate extended states appear_^ach with a localization exponent 
of « 2.3. Allowing for a finite Zeeman term leads to level attraction at high fieldstj. In fact, electron spin resonance 
Q (ESR) datac3 has shown that the spin splitting has an unusual nonlinear dependence on field. As shown below, this 
L| ' nonlinear dependence is consistent with localization phenomena and level attraction of extended states. 
• i—t'. A different type of system where two degenerate states in a strong fleld (e.g. spin states in the flrst Lan|ia|U^leY£l) 
' are coupled by random mixing, but the scalar (i.e. intrachannel) potential is absent, was recently studiecOO'EllLJ. 
5h This corresponds to, a spin-orbit coupling which is dominant relative to the scalar potential scattering. For white noise 
. . . inter-band mixingpiJ it was found that an extended-state at the original Landau level energy exists, suggesting a new 
universality class. For smooth inter-band disordered, it was suggested that there are two separate extended states as 
in the spin-split case in addition to a third extended state in between, at the original Landau level. 

In the present work we study, these and other aspects, using various extensions of the network model introduced by 
Chalker and Coddington (CC)eI. In the network model electrons move along unidirectional links which form a closed 
loop in analogy with semiclassical motion on contours of constant potential. Scattering between links is allowed at 
nodes, in analogy with tunneling through saddle point potentials in the semiclassical model. The assumption that 
each link carries current only in one direction implies that the wavepackets are sufficiently localized in the transverse 



direction, i.e. the magnetic length is small in comparison with the spacing of nodes or with the correlation length of 
the potential fluctuations. 

Our paper is organized as follows. In Sec. || we describe the network model for onCj-apd two channels and relate 
parameters of the system to the transfer matrix. In Sec. Ill , we expand our earlier workll3 and study the two-channel 
network model corresponding to two coupled Landau bands. We evaluate the energies of extended states as function 



of A, the bare energy separation of extended states in the absence of level mixing. In Sec. |V| we investigate the 
two-level system with random mixing while the scalar potential scattering is periodic (i.e. non-random) and consider 
different types of symmetries for random mixing. Our results are summarized in Sec. ^ where the ESR data is also 
discussed. In Appendix A, using the one-channel network model, we test a possibility of sublocalization behavior of 
the wavefunction, based on a possible fractal behavior of the nodes. In Appendix B the tunneling amplitude between 
Landau levels is estimated. 



II. THE NETWORK MODEL 

Consider first a one-channel CC networki, which has directed links and scattering at nodes (Fig. |l|a). Propagation 
along links yields a random phase 0, thus links are presented by diagonal matrices with diagonal elements in the form 
exp(i(^). Transfer matrix for one node relates a pair of incoming and outgoing links on the left (links 2,4 in Fig. |l]b) 
to a pair of links on the right (links 1,3); it has the form 
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(1) 



In order for the system to be invariant, on average, under 90° rotation the next neighbor node is obtained by rotating 
Fig. 03 by 90° and writing the states on links 4, 1 in terms of those on links 2, 3. The transfer matrix has then the 
same form as in Eq. with a parameter 0' replacing 9 whereO sinh^' = l/sinh0 and the phases ai, Q^a, ct4 are 
replaced by —a^, ai,a4 + tt, —a-i (Fig. |l]c). We therefore describe scattering at the nodes indicated in Fig. by circles 
with transfer matrix T(6') and at the nodes indicated by boxes with T(6''). 

The phases oli can be absorbed into phases describing propagation on links. Since a shift by a common phase in all 
links of a given column does not affect the Lyapunov exponents of the network, we can choose the phases on the links 
as ±a, ±/3 (Fig. |l|a) where a = \a\ — ^a^ — a2, l3 = ^ai — + a^. Note that the sum of left links equals those of 
right links +it. For a random potential (Section |l|) these link phases are considered as random while for a periodic 
potential (Section |v|) we choose specific sets. The appropriate procedure is then to multiply transfer matrices for 
links and nodes alternately and derive Lyapunov exponents for strips of width up to 64 and of length of typically 
60, 000 units (Section III) or 240, 000 (Section IV). At these lengths our error for the localization length (i. e. for the 
inverse of the smallest Lyapunov exponent) is ^ 0.5%. 

In the following we relate the node parameter 9 to the electron energy by using known results for scattering from 
a saddle point potentialEj. The transmission probability T of an electron with energy E through a saddle-point 
potential Vsp{x,y) — —UxX^ + Uyy^ + Vq in a perpendicular magnetic field B is given by 
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1 + exp(— 7re) 



(2) 



where 



e={E~{n + -)E2-Vo)/E^, 
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and 
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where Vt — {juj-^ + {Uy — Ux)l°ni)^^'^ , 7 = (Uy + C/a;)/4mfi, m is the electron mass and Wc = eB/mc is the cyclotron 
frequency. The oscillator frequency E2 is 
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Prom Eq. (||) the ratio of reflection and transmission coefficients is exp(— 7re). On the other hand the transfer matrix 
Eq. determines this ratio as (sinh0)^. Therefore, the relation between node parameter 6 in CC model and the 
electron energy is 

e = --ln(sinh6l) (6) 

TT 

We note that unlike discrete Landau levels, the saddle point potential allows for a continuous energy E for each 
discrete state n. Furthermore, for ^ 7 i?i ^ 0, i?2 — > t^c and E hujc{n + 1/2) + Vq corresponds to discrete 
Landau levels. In the opposite limit of Wc ^ 7 we have Ei {Ux/^my^^, E2 — > (2C/j,/m)^/^ and the integer n 
corresponds to a quantum number of the harmonic potential llyU^- 

In the case of one channel per link a single extended stateO is at e = (i.e. 9 — 0.8814), corresponding to the 
center of any band n, E ~ {n + \/2)E2 + Vq. Numerical studies of this system with width M and periodic boundary 
conditions confirm the one parameter scaling hypothesis, i.e. the localization length S^m is given by a scaling function 
/, where 

% = /f^) (7) 



M ' \£,{E)J 

with ||(-E) ~ l^'l " and — 2.5 ± 0.5. This res«]|t-|is in good agreement with experimental data for spin resolved 
levelaJ, numerical simulations using other modelst^E^ and semiclassical derivationP-E^ that predicts v = 7/3. We have 
repeated the one-channel calculation (Appendix A) with a particular emphasis on the possibility of sublocalization; 
we find that localized states decay as a regular exponential. 

A two channel network, i.e. two channels per link, is characterized by parameters 9 and 9 — A9 which determine the 
tunneling amplitude at the node for each channel with A9 related to the relatiiie energy interval between the bands 
of the two channels. The transfer matrix is parameterized in the following waylij: 

(8) 
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The transfer matrix at a node is composed of the blocks: 

f-^ _ ( cosh6' \ V. — ( ^'^^^ 1 ('Q^ 

\ cosh(6' - M) \ sinh(0 ^ M) ' " 

Defining an energy parameter A such that e — A = — (2/7r) ln(sinh(6' — /S.9) we obtain 



/ ^1 + exp[~7r(e - A)] \ / exp[-7r(e - A)/2] 

' ^1 + exp[-^e] )' \ exp[-7re/2] 



(10) 



Propagation along links yields random phases (pi,! = 1 — A and also allows mixing between two different channels. It 
is described by block U: 



u^lCi.) ^) (CI). (") 



e^P^ ) \ X \/l-a;2 ) \ e 

where is the mixing probability between different levels. 

We note that the maximal number of independent parameters in a U{2) matrix presented in Eq. ( pT| ) is four. These 
phases can be chosen as (5 = l/2(0i (t)2 + <t>?i + ?!'4),<^i = ^ l/2((/'2 + (t>i),52 = —1/2(02 — <t>?,) so that 



,5 ( e'^Wl-x^ -e'^^x 



e 



U = e-( ^ ). (12) 



The phase 5 corresponds to a scalar potential, it can be either random (Section |||) or fixed for a periodic potential 
(Section p\1). If 5 — 27rxinteger the U matrix changes its symmetry group from U (2) to a unitary unimodular SU{2). 
In Sectionpvl we study the effect of various subgroups of U{2) on the critical properties of the system. 
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III. NON-DEGENERATE LEVELS WITH RANDOM POTENTIAL AND RANDOM MIXING 



We study here a system of two Landau levels in the presence of a smooth random potential, so that for proper 
description we use the most general form of Eq. ( |l2[ ) with four random variables. 

In the absence of level mixing we know from results of the one channel systemO that e = corresponds to extended 
states. This defines "bare" extended states at E^x = E2{n + 1/2) + Vq and the bare energy splitting of the n=0,l 
states is then £2- Note that the splitting E2 which is magnetic field dominated at uJc{m/Uy^'^ > 1 remains finite 
as Wc ^ and is potential dominated at ujc{m/uy/^ < 1 (here Ux = Uy = U). The latter region is acceptable for 
a network model if the correlation length of the potential fluctuation is long compared with the magnetic length, so 
that locally the saddle point potential determines a finite splitting. The splitting ^ ^ E2I E\ is therefore bounded by 
A > 2 at ^ 0. 

The application to rt = 0, 1 Landau levels assumes that mixing of states with rt's differing by An = 2 is much smaller 
than for those with An = 1. Mixing or transition rates can be evaluated for a potential of the form (l/2)[/(a;^ + y^) — 
Uif^ / \ where A is a measure of th e co rrelation length of the random potential. For states near the local minimum we 
find that mixing is given by Eq. (B6) (see Appendix B). We therefore expect our results for the n = state, within 
the two channel model, to be valid down to rnw^/C/ « l/ln(A/^) ^ 1 (assuming U k, U) so that the range of valid 
mio\jU values is extended down to lower limits for longer range potential fluctuations. The results for the n — \ state 
are, however, not directly relevant. 

We consider below also the case of a spin split single level n = for which Vb = ±(l/2)(7*/iB, where g* is the 
electron g factor and /is is the Bohr magneton. The bare extended states correspond then to Ee,x = \E2 ± \g* \ib so 
that A = g* [IB I El. 

The system (see Eqs. (P,p^) is, on average, invariant under 90° rotation if at the next neighbor node the transmission 
and reflection (of each channel) are interchanged, i.e. e — > — e and e — A ^ — e + A. The system is then symmetric 
under e —e, A —A and the extended states at e, {i = 1,2) become ei(— A) — —ei{A). A further translation of 
energies by A returns the system to itself except for 1 2 interchange, e2(— A) + A = ei(A), hence ei(A) + e2(A) — A. 
Therefore results for the two energies £1^2 are constrained by the condition ei + £2 = A. 

We proceed further by calculating normalized localization length ^a//M; for systems with M = 16,32 (size of T 
matrix is M x M and Mi = AI/2 is the number of links across the strip) and different values of e and A. For A =p-p 
we used also M = 64 which affected critical energies by 3%, which are also within 5% from the result of Wang et al.Ej. 
The coupling x is chosen to be uniformly distributed in the interval [0, 1]; we checked that other distributions in x 
lead to similar results. Finite size scaling is then used for fitting our data onto a single curve /(Coo(e)/M;) extracting 
values of localization length ^00 for the infinite 2D system. Finally we look for the critical energies Ci and critical 
exponent v by requiring ^ k ~ £i| The raw data for one particular A (see Fig. ^) represents the characteristic 
features of the system. One can see that the values of /Mi for any two energy values whose sum equals A coincide 
as we expect from the symmetry condition. Another important feature is that for A > 1 we have two pronounced 
maxima of £,m/Mi, which we expect to be near the critical energies; this is supported by the scaling procedure. We 
find that approaching the critical energies from outside yields z/ = 2.5 ± 0.5. The states between the two critical 
energies seem to be localized, i.e. £,m/Mi decreases with M, although the decrease is rather slow. 

The critical values ei^2 are presented as a function of A in Pig-|3[ We cannot calculate ^m/M for A > 3.5 because of 
round-off errors. The A ^ 00 case can be solved analyticallyc3 since then near one extended state the other channel is 
very far from tunneling and its trajectory is a closed loop between the nodes. Eliminating the closed loop variables and 
assuming that it mixes with links of only one node leads to extended states at ei(A) ^y^/A and £2 (A) — A ^ — 1/A. 
The rather flat behavior of ei(A) up to A < 3.5 implies that the asymptotic behaviorE£l sets in at a higher A. 

The most remarkable aspect of the data is the crossing of ei(A) with the bare extended states at A « 0.5. Usually 
one expects that mixing affects mainly extended states leading to level repulsion of these states. This expectation 
leads to ei(A) approaching zero from below and £2 (A) approaching A from above, as implied by Fig. 2 of Ref. |l^. 
Contrary to this expectation we find that above A « 0.5 there is level attraction. 

Level attraction is also supported by the following obseipatipn: Parameterize the upper branch of critical e by 6, 
i.e. e = -(2/7r)lnsinh6'; at A = e = 0.23, i.e. 9 = 0.65BI13, and e increases with A so that 6* -> 0. In the level 
repulsion scenario e > A so that the parameter A9 in Eq. (^) varies in the range > A9 > —0.88. We could however 
define the model with parameters 9, A9 (instead of e. A) with A9 unbounded which implies that e — A must change 
sign. 

We present now results for the energies of extended states as a function of magnetic field by relating A to B via 
Eqs. (5,6). As noted above E2 is always finite so that A > 2 (A — + 2 for i? and Ux = Uy) so that we are 
always in the level attraction regime. The results are shown in the inset of Fig. || with the O symbol (assuming 
for simplicity Ux — Uy = U) and the full line is the lower bare extended state energy. Our data shows a minimum 
at a;c(2m/L/)^/^ ~ 0.5 in the lower state, consistent with floating, and is a result of level attraction due to Landau 
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level mixing. Allowing for mixing with the n = 2 Landau level may cause floating of the n = 1 state as well, but 
will have a small effect on the floating of the n = state, as discussed above. For N Landau levels our symmetry 
argument shows that the extended state energies come in pairs whose average is the same as for the bare states, i. e. 
Ci + EN+i-i = {N — 1) A with i = 1, . . . , N. Hence we expect that the energies of the lower half states increase at low 
fields consistent with the floating scenario and with the experimental datall3. 

The assumption of full mixing, i.e. x E [0, 1] in Eq. ( [1 1|) , is not valid for strong magnetic fields where tunneling 
between Landau levels should be suppressed. We model this situation in a reduced mixing model (see Appendix B), 
where the parameter x in Eq. (12) is now chosen randomly in the reduced range of [0, exp(— tow^/J7)]. The results 
are shown in the inset of Fig. 0, with the x symbol. The minimum in the lower level is now more pronounced and is 
at a higher field, u;c{2m/Uy^'^^ 1. 

Finally, we present application of our data to extended state energies for a Zeeman spin splitting of a Landau band 
where A = g* hbB / Ei. Since A can cross the value « 0.5 we predict that spin splitting is larger than expected for 
small fields (level repulsion) and is-amaller than expected for large fields (level attraction). In fact, we claim that level 
attraction accounts for ESR datacSI which shows a nonlinear dependence of the spin splitting on B. In particular the 
data for iV = in the range 8 — 14T is visibly nonlinear in the range 8 — 12T. 

To demonstrate the_effect we consider g* jisB /lUc = g*m/2mo = 0.018 (e.g. g* — 0.51 with an effective to free 
electron mass ratio oiEa m/mo — 0.07) and replot the data of Fig. ^ for the splitting AE = Ei (e2 — ci) in Fig. ^. The 
experimental data at B > 8T fits our results if {U/my/'^ = SlOGHz is chosen (see section V for an independent data 
leading to this value); the field scale is then ujc{U/m,y/^ « 1.3i?(T). The deviation from a straight line measured at 
8— 12T cannot be accounted for by band structure calculations (see section V). The data, therefore, provides a strong 
support for level attraction at these fields. 



IV. DEGENERATE LEVELS WITH RANDOM MIXING 



We consider in this Section a single Landau level with two degenerate spin states (i.e. no Zeeman term) where the 
only randomness comes from rpixing of the spin states, e.g. by spin-orbit scattering. A model of this type was studied 
by Hikami, Shirai and WegneiO with the Hamiltonian 

H=^ip + Af + hir)-a. (13) 



J22 



where h(r) is a random field in x, y directions and projection to states of the lowest Landau level is understot 
model has an extended state at the original Landau levelEil (e — 0) and possibly additional two extended statei 
symmetric around e = 0. 

We wish to study the Hamiltonian (Eq. ([T^)) by a network model. We therefore replace the continuum by a periodic 
potential so that the nodes of the network are the periodic set of saddle points. In general, a magnetic flux through 
a unit cell leads to Aharonov-Bohm phases on the links which increase linearly from, say, left to right. This however 
reduces the symmetry around e = as checked by our simulations; this may be related to the formal loss of invariance 
under 90° rotations. We therefore choose the unit cell to have an integer number of flux quanta. 

As shown in Figs, la-lb the phases on the links can be represented by two phases a, (3. In order to choose relevant 
values of a, (3 we consider the pure system without spin mixing and try to make its spectra similar to that of the 
p^/2m term in Eq. (13), i.e. a constant, with a total number of states close to that of a Landau level. Extended 
states in the pure system can be found by applying transfer matrices across a unit cell (i.e. two nearby nodes in Fig. 
1) which by Bloch's theorem lead to multiplication by exjp(iq) along the strip or exp(ifc) in the direction across the 
strip. This procedure yields the dispersion relation 

sing + sin(fc - a + /?) 

cosh(7re/2) = \^ . / — —. (14) 

^ ' ' 2sin(a-H/3) ^ ' 

This relation shows that there is a maximal energy emax in the band where cosh(7re/2) = l/sin(a-|-/3). The density 
of states is linear at e — s- and saturates at tmax- The total number of states in the band approaches that of a Landau 
level for a + (3 ^ Q. However, at a -I- = the spectra Eq. (14) is singular, so that it seems that a small but finite 
value of a -F /3 is needed. 

The spin term in Eq. (^3|) is taken to mix spin states on links. Since all links are unidirectional the transfer matrix is 
equivalent to an evolution operator which is a rotation in spin space exp[{i/fi) J h{r)-adt]. Since successive rotations 
about the x,y axes produce a rotation about the z axis h should be a 3D vector. The effect of these rotations 
is equivalent to an SU (2) transfer matrix (Eq. (03) with S = 0) and we can interpret all independent phases as 
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corresponding angles of the pseudofield. Below we consider such transfer matrices with different sets of independent 
parameters for the case of two degenerate levels. 

Consider first the case when T is a real matrix, i.e. the pseudofield is only in the y-direction and T has then 5*0(2) 
symmetry, 

U^f^os^ -sin</.\ 

\ sm0 COS0 / ^ ' 

The results for the case with the Haar measure, i.e. 0e [0, 27r] with cos<f> uniformly distributed, are shown in Fig. ||. 
The data shows a single extended state at e = 0. 

We know that a random U{2) transfer matrices lead to two extended states. It is therefore interesting to modify 
the 5*0(2) system by a random potential leading to a U{1) x 50(2) system 

^^^,,/cos^ -sin0\ (^g^ 

\ SllKp COS0 / ^ ' 

The data maintains a single peak at e = in this case as well. Fitting the data of those two cases by a smooth 
function / so that S,m/Mi = /{eM^^'^) yields v « 2.2 for the critical exponent (Fig. ||). Actually, one should expect in 
those two cases the behavior similar to one-channel model due to the fact that matrices U commute with each other. 
Our results support this expectation. 

The next in the hierarchy of symmetries is the SU{2) group (Eq. ( p^ with 5 = 0) with three independent phases. 
(A choice of 2 independent phases is not closed under successive transfer operations). 

The most general case of a network model corresponding to the Hamiltonian Eq. ( |l3| ) is that of SU(2) matrices on 
links and a complex transfer matrix (Eq. (^) at nodes. The T matrix on links is then a U{2) matrix, however only 
its 5*1/(2) phases are random while the phase S in Eq. ( p^ ) is regular, having the values ±a,±/3 periodically. 

Before presenting numerical results we discuss the symmetries for which ^oo(e) is invariant. Note first that a a + n 
OT P /J + TT is a symmetry, since one can shift even or odd columns by a constant phase. Shifting by tt, results in a + 7r 
and — a + 27r (the latter is equivalent to — {a + n)+n) or to /3 + tt and — /3 + tt (the latter is now equivalent to —(3 — n) . 
Consider next the up down reflection symmetry which yields (after shifts by tt) the equivalence {a, (3) — > (—a, — /3). 
The right/left reflection symmetry leads to (a, /?) — > {P+tt/2, a— 7r/2). Indeed, looking from right to left, the links with 
phases (/3, —(3) come first (actually with opposite signs since propagation is to the left, but by the previous symmetry 
signs can be changed) while (a, —a + tt) come second; shifting by ±7r this yields the stated equivalence. The next 
symmetry is a property of the infinite 2D system: if we rotate our system by 90° and then shift the phases properly 
we get (a,/3) ((a + /3 + 7r)/2, (a + /3 - 7r)/2), but due to the right/left symmetry this is just {{a + f3)/2, (a + /3)/2). 
If this symmetry holds for our finite strips then it is sufficient to consider only the diagonals in the (a, /3) plane. Note 
that all these symmetries hold for the spectra of the pure system, Eq. (14). Finally, we consider a symmetry which 
is due to the random SU (2) phases. A shift of the link phase 6 in Eq. ( [12|) by tt is equivalent to a shift of n in the 
SU{2) phases i5i, 62', since the latter are random the result is invariant. Thus in Fig. |^ the phase —a + tt can be 
replaced by —a and a column shift by 7r/2 yields (a, f3) ^ {a + tt/2, /3) or similarly (a, (3) {a,P + tt/2). Thus with 
rotation symmetry it is sufficient to consider phases {a, a) in the range < a < 7r/4. 

We checked all these symmetries numerically. It turned out that all symmetries hold, however, the rotational 
symmetry requires larger M and higher number of iterations than the U(2) case, in particular for large a, (3. This 
fact has probably the following explanation: the periodic phases a, (3 introduce a new irrelevant length scale in the 
system. In order to go beyond this scale to find the symmetries one needs to investigate much longer systems (at least 
300,000 iterations in comparison with 60,000 in the U{2) case). The rotational symmetry was found for the largest 
system size M = 128 for some particular phases. It is expected to hold for M 00. 

The results for (a = -0.175,/? = 0.075) and {a = -0.005, /3 = -0.005) are shown in Figs. 0,| respectively The 
data exhibit single peaks near e = 0. The peculiar property of the data is that in the range between e « —0.1 and 
e = the renormalized localization length ^le/S < C32/I6 < ^64 /32 which indicates a band of extended states, however 
the data for M = 128 shows ^64/32 > ^i28 /64 which determines these to be localized states. This is just another 
manifestation of the irrelevant length scale mentioned above. 

We consider finally the case a ~ /3 = in Figs. |, Although the pure case is singular, the disordered system 
is well behaved, the data converges faster and there is no irrelevant length scale, i.e. the localized behavior of e 7^ 
states sets in already at M = 32. We therefore consider a = f3 = a,s the generic case for SU{2). 

There is clearly a single extended state at e = 0, or more precisely at least all states with |e| > 0.03 are localized. 
Recall, that in the U{2) case with A = the extended states are at e « ±0.2 well separated from e = (see Refs. 
|l3|-|l5|). Hence if there are extended states in the SU{2) case near e fhO, they are extremely close to e = and 
therefore are not related to those of the U{2) case as has been proposedE3. 
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The critical exponent in Fig. is « 1.1. Thus the 5/7(2) case is a new universality class. This is consistent with 
the occurrence of singular density of states and new value of axx found by Hikami et al.Ell 

Finally, we consider an interpolation between U{2) and SU{2) by allowing the link phase {S in Eq. 12) to be random 
in a restricted range of p[— 7r,7r]. On long scales the randomness may accumulate, allowing an SU{2) behavior only 
at p = 0. However, it is known that a nietal-insulator transition (in the absence of a magnetic field) occurs at a finite 
ratio of spin-orbit to scalar randomnessE£l. Thus, it may be possible that SU{2) behavior sets in at a finite p. Fig. [ll| 
shows our data with p — 0.3; the best fit for scaling form-jpelds extended states at e = ±0.16 with v = 2.5. Since the 
extended states of the U{2) case (p = 1) are at e « ±0.2Ej"I1j we conclude that the two extended states approach each 
other as p is reduced until below p « 0.2 they merge and only e = correspond to an extended state (or states). The 
range between p « 0.2 and p « 0.05 is difficult to analyze since the localized nature of the states near e = sets in 
only at M = 128; apparently there is an irrelevant length scale, similar to the one we had above. Below p k, 0.05 the 
behavior is close to that of the SU(2) {i.e.p — 0) case, i.e. Figs. p|JTo|. Thus we have a J7(2) to SU{2) phase transition 
at a finite p, p ~ 0.2. 

We summarize our results in Table I. The table shows the different symmetries involved with their critical properties. 

V. DISCUSSION 

We have studied two types of systems: (i) Non-degenerate states with random mixing and random scalar potential, 
and (ii) degenerate states with randomness only in the mixing terms. The system with non-degenerate states is 
relevant to two types of experimentally studied cases. The first case is where non-degeneracy is represented by 
Landau level splitting. W e fin d that the lower Landau level has a minimum as function of magnetic field, consistent 
with the floating scenaricOO. This result-accounts for a transition from a Hall liquid to an insulator at both high 
and low fields, as observed experimentallyllj. 

The second case is where non-degeneracy is represented by the Zeeman spin splitting. In this case the bare splitting 
A (which is g* hbB / Ei in the simplest case) can cross the value « 0.5 where level repulsion crosses into level attraction. 
We predict therefore that the spin splitting of the extended states is larger (smalleij)-.than the bare Zeeman splitting 
at low (high) fields. To estimate the field corresponding to A w 0.5 we use dataEZi on temperature scaling of the 
conductance peak in GaAs / Al^Gai^xAs which shows that the low temperature scaling of the N = Q [ peak neax 
B = 6T is modified at a crossover temperature of w Q.iK . Analysis of tunneling through saddle point barriersEj 
shows that the crossover temperature is Ti — hU / (pjrmujc) where Ux = Uy) is the barrier curvature. Since 

Ti < ~ 12QK at B ^6T (using an effectiV|&,massE!l of m/mo ~ 0.07) Eq. (4) yields Ei « TiU/mLOc = 2itTi « 2K 
and (fiU/mY/'^ ~ 15K. Using the bulk valueEBl g* — 0.44 we obtain Awl, i.e. level attraction is predicted above 
B^4T (where A « 0.5). 

Spin splitting of extended states can be directly probed by ESR; the extended electromagnetic wave couples domi- 
nantly the extended states. Experimental dataw on Landau levels iV = 0, 1, 2 ahjpws a spin splitting which is nonlinear 
in B and was fitted by a quadratic polynomial. Band structure calcuiationsESo show that at large fields the 
system has a spin splitting of the form AEq ~ g* fisB + i^n, assuming23 z^jv ^ ^c', (in the more recent calculationEa 
the linear form is valid for N ^ 0, B > QT). 

The nonlinearity is best seen in the N. ^ data which spans B = 8 — 14T and is clearly nonlinear in the range 
8 — 12T. The band structure calculationstjo predict linear behavior, at least at _B > 6T. We propose then that the 
deviation from linearity is due to localization effects, i.e. level attraction. Npte also that if indeed the ESR transition 
probes only extended states, it accounts for the remarkably low linewidthEHl of 50mT. _ 

To fit the data we choose AEq = g* ^bB + va with g*m/2mQ — 0.018 and i^o = (i.e. forE3 m/mo — 0.07 we have 
g* = 0.51). The parameter A of Fig. 3 is now A = AEq/Ei, allowing the data of Fig. 3 to be replotted as the 
splitting AE = i?i(e2 — ei) in Fig. 4. The experimental data could also be fitted with a smaller g* by increasing vq, 
e.g. g* = 0.44 with vq = llGHz. The splitting in frequency units is given by (C7/m)i/2 = 310GHz while the field 
scale is set by uJciU/my/"^ = 1.3B(T). r-| 

Fig. 4 shows nonlinear behavior with a decreasing slope for B > 8T, in agreement with the experimental datac3. 
Below 8T we predict a different trend with the energy splitting crossing the bare one at « 4T (Fig. 4) and becoming 
larger than the bare splitting at _B < 4T. (Recall that the bare splitting itself, which should be determined by band 
structure calculation, may be nonlinear in B for these lower fieldaH3o). We propose then that extending the ESR 
data to lower and higher fields can serve as a unique tool for testing our predictions for level attraction and repulsion. 

We have also studied the network model with random mixing of two degenerate states, while the scalar potential 
is periodic (non-random). The random part of thp transfer matrix Eq. (12) is then an SU(2) matrix. We believe that 
this model corresponds to that of Hikami et al.tiJ Eq. (13). There are, however, some differences: a) The continuum 
translation symmetry of the free particle term in Eq. (13) is replaced by a periodic potential with extended states in 
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the band of Eq. (14). This band resembles the continuum one for small a, f3. b) The random fields h^, hy of Eq. 
(12) correspond in our case to a general SU{2) matrix, i.e. we have also a random hz component. While Eq. (13) 
including a random hz was not explicitly studied, we believe that the noncommutative Pauli matrices ax , cry generate 
az terms in the time evolution (or in perturbation theory) and the models are therefore equivalent. 

Our results for the SU(2) model show that for small but finite a, (3 a single extended state appears at e = 0, though 
the localized nature of nearby states appears only at unusually wide strips {M — 128). This indicates that the phases 
a, (3 produce a new (irrelevant) length scale which renders numerical analysis more difficult. Remarkably, for the 
choice a = P = the irrelevant length scale disappears and we find a clear single extended state at e = 0. Finite 
size scaling is well obeyed (Fig. ^0|) and an exponent i/ = 1.1 is found. Thus the SU(2) network model forms a new 
universality class. |— ■ P- ■ 

Our SU (2) data are inconsistent with Lee's arguments and with the numerical study of Minakuchi and HikamiEZl. 
Lee's argument is qualitative and applies to slowly varying random fields h(r) while in our case h(r) is uncorrelated 
between unit cells. We do not expect that this will change critical properties since correlations will vanish along long 
loops in both cases. The numerical data of Minakuchi and Hikami, who study Eq. (13), shows 3 extended states. At 
present we cannot account for this discrepancy with regards to the location of extended states; although Eq. (13) is 
defined differently from our SU{2) model, the two should be equivalent, as discussed above. 

Experimental realizations of the SU{2) model correspond to systems with dominant spin-orbit coupling. We predict 
that in a usual random system which has spin split levels, increasing the strength of spin-orbit coupling would shift the 
extended states so that they approach each other, until at some critical value of the spin orbit coupling, the energies 
of these two extended states would merge. We simulated this situation by allowing randomness in link phases to vary 
in the range p[— tt, tt] range. We find that below p k, 0.2 the system behaves like an SU{2) one, i.e. an extended state 
(or states) at the single energy of the band center. 
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APPENDIX A: TEST FOR SUBLOCALIZATION 



ased on the semiclassical picture. An electron with 
on clusters of size of the percolation length which 



The simplest argument for deriving the 7/3 exponenlu is 
energy E follows semiclassical constant potential contours 
according to percolation theory diverges as ~ jiJl"*/"^. An electron can propagate to a distance r ^ by tunneling 
events between clusters, each event reducing the wavefunction by a factor exp(— a|i^|) and a depends on potential 
barrier parameters. These tunneling events occur at saddle points in the potential which are on the infinite percolating 
cluster at = 0. Assuming that these saddle points are homogeneously distributed, i.e. their number is r/^p, we find 
that the total wavefunction decay is 

[exp(-a|£;|)]"/«^ =cxp(-r/C). (Al) 

Hence the localization length f ^ (,p/\E\ ^ \E\^'^/^. 

We wish to test the homogeneity assumption, i.e. whether the number of tunneling events is r/^p. In the theory 
of percolation the infinite cluster at the critical energy has a backbone consisting of alternating sequence of singly 
connected ("red") bonds and multiply connected ("blue") bondstHTEj. It is conceivable that tunnelieg occurs at "red^ 
bonds of the infinite cluster, i.e. bonds whose elimination disconnects the infinite cluster. AnalyticE3 and numericaS 
studies show that the number of " red" bonds L^ed between two points separated by a correlation length ^p diverges as 
E~^; therefore at the scales r = ^ E~^^^ we obtain Lbb ~ r^/'*. If we assume that the number of tunneling events is 
(r/^p)'^/^, the derivation of the localization length Eq. ( |A1[ ) is now modified to exp[(— ai?)(r/^p)^/'*] ^ exp[— (r/^)^/*]. 
The consequences are that now the localization length exponent is modified ^ ~ E~^^^, and more significantly that 
sublocalization appears with a weaker exp[— (r/^)^/^] decay. To check this possibility consider how it affects finite size 
scaling and in particular the region of localized statesLj. Suppose that wave function on an Ad x Ad square decreases 
as exp[— (AZ/^m)'']- Then at distance r > AI, where M is the width of the strip, it decays as in a ID strip, i.e. as 
exp[—{AI/£_M)'^{r/AI)] ^ exp[— t/^m], and is the value that we calculate numerically; hence, £_m — {£.m/M)^A4. 
On on the other hand, by scaling hypothesis ^m/M — f{(,/A4), which for ^ <C M yields £,m/M = £,/Ad and therefore 
£,m/M = {(^/AiY . Therefore an asymptotic tangent to the fitting curve of log(^M/-W) versus log(M/^) should have a 
slope —fi (—3/4 if the idea of "red" bonds is valid). We perform optimization procedure using our results as well as 
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data from CCEj. The fitting curve (Fig. |T^) and its tangent in the range ^ <C M show conclusively that /? can not be 
distinguished from unity and therefore the localized wave functions decay exponentially. 



APPENDIX B: MIXING OF LANDAU LEVELS 

We wish to estimate the mixing probability of Landau levels on the links, i.e. far from saddle points. We consider 
then the tunneling rate between Landau levels in a weakly perturbed harmonic potential. We choose the potential 
on the links as V{x,y) — {l/2)U{x'^ + y^) — {U/X)y^; A is a a measure of the correlation length of the potential. We 
perform a standard change of variables 

7r = p-— , X Y —. (Bl) 

X and Y are slow guiding center coordinates and T:x,y are operators of fast cyclotron motion. A Schrodinger equation 
for the fast variables (choosing ctIx/ eB = z as a coordinate and treating X and Y adiabatically) has a form 




2,UY\ ( c7r/m - (3i7/Am)y2 \ tj ^ 
\m j y oj^ + jj/m- {6U/Xm)Y J 

xV„(z;X,y) = ^„(X,r)^„(z;X,r). (B2) 



Eigenvalues of Eq. (B2) can be interpretedE3 as local nonequidistant Landau levels En{X,Y) depending on the 



classical guiding center coordinates. 

We estimate the tunneling rate between neighbor Landau levels using Dykhne's formulatllH 

P„,„+i ~ exp[-(2/?i)Im / \e,,+i ~ E,,)dt], (B3) 

where En+i > En for all real times, and tc is a point in the complex time plane where they cross. We solve Eq. ( ^^ ) 
keeping X and Y as fixed parameters and obtain to the first order in 1/A 

£;„+! ~En^nu;,{l + -^~ 3UY/X^ 

Note that we are considering quasi bound states near the local minimum of V{x, y), i.e. Y <^ X. We solve equations 
of motion for 1" in a harmonic potential, substitute this solution into Eq. ( |B4| ) and find 

TTTOWc , .muc , 2X{mujl + UY 

tc ^ ^zbi^^ln — ^ ^ —, (B5) 

2 U U imUulA 

where A is the amplitude of oscillations in the F-direction. Note that UA is a measure of the center coordinate 
energy which corresponds to the width of a Landau band; since this width is less than or on the order of TloJc we have 



A < \ TiWc/U — idrnuol/il where ^ is the cyclotron length. We find crossings at conjugate points as one expects for 



real HamiltoniansH. Evaluating then Eq. (B3) we find that the main contribution to the tunneling (in this adiabatic 



approximation) is independent on n and proportional to 

/ _, ,2 , jj 2X{mujl + Uf 



'„,„+! ~ exp 1^ In ^ ^ ' . (B6) 



For transitions from n to n + An the term in the exponent is multiplied by An. 
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FIG. 1. The network model. Arrows indicate direction of current flow. The nodes marked by o and □ are related by a 90° 
rotation, (a) The strip has Mi links with periodic boundary conditions {Mi = 6 in the figure). In the two channel case each 
link represents two channels, i.e. M = 2M; is the number of transverse channels. The thick line is a reference for measuring 
area in Section 4. (b) and (c) Phase relations for two neighbor nodes. 
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FIG. 2. Renormalized localization length (^m/M as a function of e for A = 2.2. + symbols correspond to M = 32 and 
□ to M = 64 system widths. Arrows point to the location of extended states in the absence of level mixing. The energies of 
extended states (near the peaks of (,m/Mi) are closer than the arrow positions, demonstrating level attraction. 
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FIG. 3. Critical values e as functions of A. Full lines arc the bare extended states e = and e = A are shown. Inset: Energy 
of the lowest extended state as a function of magnetic field with full level mixing (O) and with reduced mixing (+). Full line 
is the lowest bare extended state. Here E' = E{2m/Uy^^, oj'^ = Uc{2m/Uy^^. 
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FIG. 5. Renormalized localization length (^m /Mi as function of energy e for a system with SO{2) mixing between the levels. 
O symbols correspond to M = 16, 4- to M = 32 and □ to M = 64 system widths. 
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FIG. 6. Fit of raw data from Fig. g for a system with 50(2) and for a system with f/(l) x 5*0(2) mixing between the levels: 
^m/Mi versus eMl^" with = 2.2. O symbols correspond to M = 16, + to M = 32 and □ to M = 64 system widths for 50(2) 
case; x symbols correspond to M = 16, A to M — 32 and ★ to M — 64 system widths for (7(1) x 50(2) case. 
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FIG. 7. Renormalized localization length £,m /Mi as function of energy e for a system with SU{2) matrices on links and and 
a complex transfer matrix at nodes with small (a — — 0.175, /3 = 0.075) phases after 480,000 iterations. O symbols correspond 
to M = 16, -I- to M = 32, □ to M = 64 and x to A/ = 128 system widths. 
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FIG. 8. Renormalized localization length /Mi as function of energy e for a system with SU (2) matrices on links and and 
a complex transfer matrix at nodes with small (a = /3 = —0.005) phases after 480,000 iterations. O symbols correspond to 
M = 16, + to M = 32, □ to M = 64 and x to M = 128 system widths. 




FIG. 9. Renormalized localization length (,m/Mi as function of energy e for a system with SU{2) mixing between the levels. 
O symbols correspond to M = 16, -|- to M = 32 and □ to M = 64 system widths. 
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FIG. 10. Fit of a raw data from a Fig. |^ for a system with SU{2) mixing between the levels: ^m/Mi versus eM^" with 
u = 1.1. 
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FIG. 11. Renormalized localization length £,m/Mi as function of energy e for a system with reduced (7(2) mixing between 
the levels (p — 0.3). O symbols correspond to M = 16, + to AI = 32 and □ to AI = 64 system widths. A fit to a scaling form 
yields two extended states ai e — ±0.16. with v — 2.5 
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FIG. 12. A fit of the data for one-channel model: £,m/M versus M/^, and a tangential straight line with slope equal to unity. 

TABLE I. Symmetry classes of random phases 

Group link phases a, /? # of energies with extended states v 

SO{2) arbitrary 1 2.2 

SU{2) 1 1.1 

U{2) random 2 2.5 



18 



